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1 Introduction 

The main purpose of this note is to give a simphfied version of Lieb's proof of the flux phase 
conjecture in which at the same time allows for some straightforward generalizations. 
Those readers who are mainly interested in the basic argument, rather than in learning 
about the more general description of it, are advised to think about a finite regular square 
lattice on a cylinder while reading this and the next section. Once the argument is properly 
understood the generalizations become straightforward. 

The physical context where the first conjectures appeared |@] is reviewed in Q. For a 
history and a more general formulation of the problem we refer to the first mathematical 
studies on the subject by Lieb and Lieb and Loss 

Consider a system of spinless fermions (adding spin poses no problems) on a finite set of 
sites A, at half-filling, and with a Hamiltonian of the form 

H= Y: txyclcy + H^^^ (1.1) 

x,?/GA 

Here, t^y is a hermitian matrix. We will explain later what kind of interactions H^^^ are 
allowed (see Section ^). For now, let us just say that the usual on-site Hubbard interaction 
is among them (i.e. for spin 1/2 fermions one takes H^^^^ = UJ2xeAi^x,i — 1/2) (^x.i — 1/2) ), 
and that only gauge invariant interactions will be considered. 

Let F be the graph with set of vertices A and the set of edges {< x,y >\ t^y ^ 0}. 
A circuit in the graph F is a finite sequence 7 = (xi, . . . , Xk) of distinct vertices, with the 
property that < Xi, Xj+i >, for i = 1, . . . , — 1, and < Xk,xi >, are all edges in the graph. 
By representing the circuit as an ordered sequence we have implicitly given it one of the two 
possible orientations (for k > 2). 

The ground state energy of ( |1.1| ) depends on the, in general complex, parameters t^y only 
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through their modulus \txy\, and the flux variables for circuits 7, which are defined as 
follows: 

n 

*7 = IZv'x^xi+i, mod27r (1.2) 

i=l 

where 7 = {xi, . . . ,Xn), and t^y = exp{iipxy)\txy\- This follows from Lemma 2.1] where 
it was proved that there is a unitary transformation relating any two Hamiltonians with 



phases {(fxy \< xy >E T} that satisfy ( |1.2| ) with the same fluxes $-y for all 7. This unitary 
transformation is of the form 

cl^e^'-cl c^^e-^'^c, (1.3) 
and is called a gauge transformation. We will often write {^Pxy} instead of {(pxy \< xy >G F}. 

The flux phase problem can now be formulated as follows: for fixed values of the moduli 
\txy\, find the phases ipxy (or, equivalently, the fluxes $^) for which the ground state energy 
of the Hamiltonian ( |1 . 1| ) attains its minimal value. We cannot solve this problem in general. 



In fact, we do not expect that there is a simple solution in general. We are looking for a 
solution in terms of a basic set of circuits C (e.g., the plaquettes of the square lattice). The 
set C should be not too large and consist only of "simple" circuits, so that the solution (= 
the values of the fluxes through the circuits in C) can be easily described. On the other hand 
C should contain enough circuits so that their flux uniquely determines $^ for all circuits 7. 

Definition 1.1 A set C of circuits in a graph T is called a basic set of circuits if for any 
two configurations of phases {fxy} o^nd {f'^y} that produce the same fluxes | 7 G C}, 
there exists a gauge transformation relating {^fxy} and {^p'xy}, i-e. ip'^y = (fxy + Oy — Ox, for 
some real 6x,x G A. 



Lieb and Loss showed in P, Lemma 2.1] that the set of all circuits of a graph satisfies 
Definition LI. Often, it is more convenient to work with a rather small subset of the set of 
all circuits. For examples of good choices of the set C we refer to Section ^ 
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The class of models that we treat in this paper is described by the following two as- 
sumptions on the graph T together with the configuration of \txy\ 's associated with the 
bonds: 

Al. All circuits 7 = {xi, . . . , Xn) in F are of even length, i.e. n — 2k. This is equivalent to 
requiring that the graph F is bipartite, but we will not use explicitly a decomposition 
into two sublattices. 

A2. There is a basic set of circuits C such that for each 7 e C there is an embedding of the 
graph in R^, for some D, such that there is a D — 1-dimensional reflection hyperplane 
P not containing any vertex of A, with the following properties: 

1. The whole graph F, together with the configuration of l^aij/l's, is invariant under 
reflection through P. 

2. All circuits 7 e C that are intersected by P (i.e. not all vertices are in one of the 
two halfspaces) are, up to orientation, invariant under reflection through P. In 
particular, 7 is invariant under reflection through P. 

The embedding of the graph in used to describe Assumption A2 is not essential. We 
only introduce it in order to simplify the description. 

Before we can state our main result we have to say what we mean by "flux configuration" 
and "canonical flux conflguration" . 

Definition 1.2 (Flux configuration) Let {$7} be a set of fluxes (i.e. real numbers mod 

27r) through all circuits of the graph. We say that {$7} is a flux configuration, if there exist 
a set of phases {<Pxy} such that (1.2) holds for all 7. 
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Figure 1: Two different embeddings of a grapli in tlie plane. Tlie solid lines 
indicate the edges. The dashed lines show a triangulation. The Lieb-Loss flux 
through the circuit (1,2,3,4,5,6,7,8) is for the first and tt for the second 
embedding. 



Definition 1.3 (Canonical flux configuration) Assume that T satisfies the assumptions 
Al-2. A flux configuration is called canonical if there is a set C of basic circuits 

satisfying A2 and such that for all 7 G C, $^ = z/ 7 has length 2 mod 4, and $^ = vr z/ 7 
has length mod 4. 

Note that it is not true, in general, that in a canonical flux configuration all circuits 
satisfy $^ = if 7 has length 2 mod 4, and $^ = tt if 7 has length mod 4. 

The arguments in Section |^ will show that, for graphs that satisfy the assumptions Al-2, 
there always exists a canonical flux configuration. 

The definition of canonical flux configuration given here is different from the one put 
forward in for planar graphs embedded in the plane. A planar graph embedded in the 
plane can be triangulated and Lieb and Loss [^] note that the number of triangles enclosed by 
a circuit is independent of the triangulation, and they define a flux configuration by putting 
a flux 7r/2 in each triangle. The resulting flux conflguration for the original graph, however. 
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depends on the embedding in the plane one starts from (see Figure 1 for an example). Our 
definition is restricted to graphs that have a basic set of circuits satisfying the assumptions 
Al and A2. They need not be planar, but, on the other hand many planar graphs do not 
have a canonical flux configuration according to the definition given here. Also we do not 
know whether there are graphs for which different choices of C lead to different canonical 
flux configurations. 

Our main result is the following theorem. 

Theorem 1.4 Under the assumptions Al and A2 we have the following: 

i) There exists a configuration of phases {v^^y} such that the corresponding configuration of 

fluxes is a canonical configuration. 

a) For the Hamiltonians ^i. j| j we have 

inf UH{{v.^y})) = Ao(i/({^g})) (1.4) 

where Xo{H) denotes the smallest eigenvalue of H, i.e., canonical flux configurations mini- 
mize the ground state energy. 

Quite generally we expect the energy minimizing flux configuration to be unique up to 
gauge transformations, but we have not studied the question of uniqueness. Non-uniqueness 
could arise in two ways. If there is more than one canonical flux configuration the minimum 
will be attained in both. The other possibility is that there is a non-canonical minimizing 
flux configuration. 



On the Qux phase conjecture 



7 



2 Proof of the Main Result 



First, we only consider non- interacting spinless fermions. The Hamiltonian is (|1 . 1|) with 
-^int ~ 0- ^ill indicate at the end of this section how spin and certain interactions can 
be included. 



Statement i) of the Theorem |1.4| , for the case of planar graphs, is a consequence of [§[ 
Lemma 2.2]. For the more general situation considered here i) will be a byproduct of the 
proof of ii). 

The main argument is an application of the Dyson-Lieb-Simon Lemma in the following 
form. 



Lemma 2.1 Let A, B,Ci, . . . ,Cn be a collection of d x d complex matrices (n could he 
infinite) with the following properties: A and B are Hermitian, and for all i, Ci is real and 
X]j Ci ® Ci is symmetric (as a d"^ x d"^ matrix). Let Xo{A, B) denote the lowest eigenvalue of 
the matrix 

T{A,B) = A(^l + 1(^ B -Y,C^®C, (2.1) 

i 

Then 

Ao(A, 5) > 1 (Ao(A, A) + Ao(5, B)) (2.2) 

where A denotes the matrix obtained from A by complex conjugation of the matrix elements. 
In particular 

Xo{A,B)>mm[\oiA,A),\o(B,B)) (2.3) 

In the formulation of this lemma in the matrices A and B are required to have real 
matrix elements. It is crucial for our application that we consider complex matrices A and 
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B. This is a straightforward extension. For a proof of Lemma p.l| in the zero-temperature 
form stated here see 0]. 



Before we can apply this lemma we have to bring the Hamiltonian into the form (2.1). 
This will be achieved in three steps each consisting of an elementary transformation. 

Given a circuit of the set C we consider an embedding of the graph in and a reflection 
plane P of the circuit (it exists by assumption). This defines a left part (L), a right part 
(R), and a set (M) of vertices which belong to edges < x,y > with x E L and y E R oi 
X E R and y E L. The three steps are: 

i) A Jordan- Wigner type transformation, 

ii) A particle-hole transformation, 

iii) A gauge transformation. 



We know from experience that one easily gets confused while performing this sequence of 
transformations. Therefore, we now spell them out in detail and indicate the purpose of 
each of them. 

Step i): We introduce new operators df defined by 

«!. = (-ir^c,, 4 = 4(-l)^^ (2.4) 

for all a; G A, and where is the total particle number in the left half of the lattice, i.e., 
= J2xeL clc^- If one considers fermions with spin, Nl has to be the total particle number 
on the left, i.e., = J2xeL,aclaCxa- In one dimension, the transformation defined in ( p^) is 
similar to the usual Jordan- Wigner transformation. Strictly speaking however, even in one 
dimension, it is different. A slightly different transformations was employed previously by 
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several authors, e.g., in |^. Note however that in the paragraph about fermions contains a 
mistake. With the transformation employed there the hopping terms on the right acquire the 
opposite sign of the hopping terms on the left, and thus the Hamiltonian is not in reflection 
positive form. 

Using the canonical anticommutation relations of the c operators, one easily finds that 
the d operators satisfy the following algebra: 



The operators df acting on Fock space (associated to A) can be identified with operators of 



acting on the tensor product space 'Hl'^'Hr, where each factor corresponds to the Fock space 
associated to the left and right parts of the lattice. In terms of the df the Hamiltonian can 
be considered as acting on Hl ® Hr and takes the form 






the form 



d* for X e L 

l0d* for X e 



H 



X/ ^xydxdy 



(2.5) 



E t^ydxdy+ E txydldy+ ^ t^ydld.^ 



(2.6) 



x^ydL x,y(iR x,y^A'I 



The third term of ( p.6| ) describes the interaction between the left and the right half of the 
lattice and is of the tensor product form as in ( p.l| ). 
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Step ii): The second step is a simple particle-hole transformation on the right half of the 
lattice, i.e., for all x e it! 

dx-^dl dl^d^ (2.7) 
while the df with x & L remain unchanged. The Hamiltonian becomes 

^ = E txydld^+ E {-t^y)dld^ (2.8) 

x,y£L x,y£R 

+ E t^vdldl+ E t^ydxdy (2-9) 

xeL,y£R x€R,y&L 

Step in): Finally we perform a gauge transformation with the purpose of making the hopping 
matrix elements across the reflection plane all negative. A transformation that achieves this 
is the following: 

dl ^ -e-'^^-ydl 
^ " (2.10) 

dy ^ -e'^^^ydy 

for sites y E R which are connected to a site x e L (i.e. given y & R there exist an a; e L 
such that txy 7^ 0). 

Therefore we have a new set of hopping matrix elements {t'^y} with the same fluxes as 
the original configuration {txy}, (because a gauge transformation does not change the fiuxes) 
and \t'^y \ — \txy\, in terms of which the Hamiltonian is 

^ = E i'.ydldy+ E i-t^y)dldy (2.11) 

x,y&L x,yGR 

- E ICl4<- E KKdy (2.12) 

xeL,yeR xeR,yeL 



We denote by $l (respectively ^r) the set of fluxes through basic circuits which are 
entirely in L (respectively R), and by the flux conflguration for basic circuits which 
have the same reflection plane P. These fluxes refer to a particular orientation of the basic 
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circuits: first orient in an arbitrary way all circuits on the left, and for the circuits on the 
right take the orientation opposite to the one obtained by reflection of the left part. For the 
ones in the middle choose an arbitrary orientation. For a fixed configuration of \txy\, the 
ground state energy depends on the phases ipxy only through the fluxes and we will denote 
this energy by Eo{^l, ^m, ^r)- 

We adopt the convention that the same set of fluxes $i when it appears as the third 
argument of Eq, assigns the flux to a circuit on the right that is associated by reflection to 
the circuit on the left. We can now state the basic lemma. 

Lemma 2.2 Assume that the configuration {\txy\} is invariant under reflections. Then 



where $y is the canonical flux configuration through the basic circuits intersecting P. 



(2.13) 



proof: 



The proof is a direct application of Lemma 2J to the Hamiltonian in the form ( |2.12| ) , while 



carefully keeping track of the flux configurations. The operator T{A, B) of ( p.!]) is given by 



x,yGL 



B = i-t'xy)dA 



x,yGR 



= E \t'xyKdl+ E KyKd, 



y 




and Ao(A, A) is the ground state energy of the Hamiltonian 



nA,A) 



E ^xydidy+ ^ t'xyd\(x)dr{y) E l^yl'^I'^y E l^yMx'^; 



y 




E ^'xydldy+ E ^'xydl(y)dr{x) - E I^^^xj/MI^ ~ E 



y 
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2n 1 




Figure 2: By assumption, when a circuit is intersected by the reflection plane, 
it is reflected into itself. Such circuits constitute the "middle part" of the 
graph. The figure also shows a circuit (1,2,3,4) in the left part of the graph 
and its reflection (r(l), r(2), r(3), r(4)) on the right. 

where r{x) denotes the reflection of the site x through P. This can be written back in the form 
of ( |2.12| ) with a new configuration of hopping matrix elements t"y which do not, in general, 
have the same fluxes as the original hoppings. We now determine the new configuration of 
fluxes in,ni,^R)- 

First, take a circuit in the middle part 7 = [xi, ...,X2n)- We label the vertices such 
that xi, . . . ,Xn € R, and Xn+i, ■ ■ ■ , X2n ^ L (see Figure 2). The edges intersected by P are 
< X2n-,xi > and < >. The corresponding term in the transformed Hamiltonian is 

2n-l ra-l 

^Xi.Xi+i'^Ii'^ajj+i + ^ (~^x.,a;,+ i)^I,^2:i+l (2-14) 
i=n+l i=l 

-|C..ML< - IC+i,.j4„+,4„ + h.c. (2.15) 

and the corresponding flux is the original one = I]i=i fxi,Xi+i = J^iZn+i f'xi,Xi+i + 
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y]"_i if' „ mod 271. After reflection it becomes 

2n-l 2n-l 



=n+l i=n+l 

-|C..ML< - IC+.,xj4„,.4„ +h.c. (2.17) 



The new flux on 7 is = Ei=„ii '/'L^.^.+i + E?=Ji(vr - v'x,,x,+,) mod 27r, which is equal to 
(n — l)7r, i.e., the canonical flux through the circuit 7. 

Next, we consider a circuit on the left 7 = (xi, ...,X2„), oriented from xi to X2„ (n is an 
integer), and its reflection r(7) = (r(xi), r(x2n)) on the right oriented from r(x2n) to r(xi) 
(See Figure 2). After the transformations i)-iii) the corresponding terms in the Hamiltonian 
are (with the convention X2n+i = xi) 

2n 



A = E4,..+i<4,+i+h.c. (2.18) 
1=1 

2n 

B = E(-<(..,.),Kx.))4(..,.)^Kx.)+h.c. (2.19) 

i=l 

and the fluxes through 7 and r(7) are respectively $^ = 1]^=! 'fxi,xi+i = S?=i V^x^.x^+i 
'^'r(7) = E?=i V5r(z,+i),r(x.) = E?=i V5r(x,+i),r(x,) " When wc apply Lemma 2.1 we have to replace 
B by 

2n 2n 

^ = E C^^I(x,)^K-.+i) + h-c. = E t'xux.+^dl,+,,)dr{x,) + h-c. (2.20) 



i=l j=l 



The new corresponding flux through r(7) is $'/(^) = Z]i=i(7i" ~ which is equal to 

— In particular if 7 on the left has the flux or vr then r(7) on the right has the same 
flux. 



One argues similarly for Xo{B, B). This ends the proof of Lemma ^l2 . 



Proof of Theorem \1.4\ By assumption the configuration {\txy \ } is invariant under reflections 
through all reflection planes of the circuits in C. The crucial property is that for each basic 
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circuit there is a reflection plane that intersects it and for which the conditions of Lemma 



are satisfled. The theorem is then proved as an apphcation of Lemma p.2| . The Lemma 
yields the existence of a conflguration of fluxes for which the ground state energy is at least 
as low, while at the same time the new flux conflguration is produced from the old one by 
either 

(<l>i, ^ L, or (<l>i, ^ $ (2.21) 

In both cases the flux in all circuits intersected by P becomes canonical. By the same 
argument as in 0, Proof of Theorem 4.2] or one can now prove that the minimum is 
attained in a canonical configuration by showing that, in an energy minimizing configuration, 
the maximum number of circuits in C with canonical fiux must be the total number of circuits 
in \C\. Let {$7} be a minimizing configuration with a given number Ncii^^y}) of circuits 
(in C) with canonical fiux, and let 70 G C be a circuit that does not have canonical fiux 
in that configuration. Let P be a refiection plane leaving 70 invariant. After refiection the 
new configurations in ( |2.21 ) both have the same minimal energy. Then writing = 



($i,$M,$ii) 

iV,($^, + iV,(-$K, $g, = 2 (iV,($i, $h) + iVc($!S) - iV,($. 

As 7o is a circuit in M that does not have canonical fiux in $m while in it does (just 
like any other circuit of C intersected by P) it is clear that Nc{^^^]) — Nc{^m) > 1. We 
conclude that at least one of the new minimizing configurations has strictly more circuits 
with canonical fiux than {^-y}. This argument is then repeated until all 7 G C have canonical 
fiux. ■ 

Remarks: 



a) Finite temperatures. Lemma 2.1 holds with Aq replaced by — logtr exp(— ). Thus 
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Lemma |2.2| holds also with the ground state energy Eq replaced by the free energy (at half 
filling), and of course its proof and the proof of theorem |1.4| is the same. 



b) Interacting systems. It is straightforward to generalize the proofs to include spin and 
some class of interactions. One can accommodate for example a Hubbard term 

xeA xeL x&R 

where = Uijix^ — \){nxi — |), n^a = c^^Cxa, o" =1; i; U is an arbitrary real number. 
Another example is a nearest neighbor repulsive potential 

= ^ + ^ hxy + ^ hxy 

x,y£A x^ydL x,y£M x,y£R 

where h^y = V{nxi + n^i — l)(?T-j/t + ^yi ~ 1) with V a positive number. Longer range 
interactions, and spin dependent forces such as a Heisenberg antiferromagnetic exchange 
term can also be included. These cases are also discussed in |jl|. Let us describe what 
happens in the transformations i)-iii). In the first step i) one has to replace ( p.4| ) by 

d = (-l)^^c , =c^ (-1)^^ 

X(T V / X(T ' X(T xa\ I 

with Ni^ = J^xeii^xi + nxi). In step ii) for x G -R, n^a — 1 — n^a- The Hubbard term 
remains unchanged, but the nearest neighbor interaction becomes 

^xy = ^ hxy — ^ hxy + ^ hxy 

x,y£A x,yGL x,y£M x,yGR 

Thus the interaction between the left and right parts of the lattice is of the form J2i Ci ® Cj 
with the correct sign because ^ > 0. The third step iii) is a gauge transformation which does 
not affect the interaction terms. Summarizing we see that we can bring the Hamiltonians 
into the form (|2.1|) . Then the proofs of Lemma and Theorem |1.4| are unchanged. 



We believe that these remarks are useful in other problems. We illustrate this by two 
examples: the t — V model and a generalized Falicov Kimball model. 
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c) Spinless t — V model. This model of spinless electrons has Hamiltonian |1.1| with the 
interaction part equal to V J2x,yeA{^x — ^/'^)iny—l/2) where the sum is over nearest neighbors 
only and V is positive. The remarks above show that on a cubic lattice, i.e., D = 3, with 
periodic boundary conditions in all directions and a flux configuration through each square 
plaquette equal to vr, and \txy\ = t for all bonds < >, it can be brought in a reflection 
positive form with respect to all reflection planes. Then it is an exercise to see that the 
methods of used for the Heisenberg model can be used also in the present situation to 
prove that, when t/V is small enough there is long range order at low temperature (3~^. 
More precisely if < ■ >a is the thermal average with periodic boundary conditions, one can 
prove (— l)l^l+ls^l < [rix — l/2){ny — 1/2) >a> c > 0, for all x and y in A, for some strictly 
positive constant c independent of A. (We note that in the present case the uniform density 
theorem ||10| applies, so in particular < Ux >a= 1/2 for all /5, t and V). In fact this result 
is true for any flux configuration and one can also add a small chemical potential term (see 



TT|],[|T^ and for recent rigorous results). Although our proof does not work in D = 2, 



the result is expected to hold also in two dimensions. 

d) One can also consider the case of spin 1/2 electrons with attractive Hubbard interaction 
and a nearest neighbor repulsion, i.e., a.t — V — U model, and it can be shown that for low 
enough temperatures, in three or more dimensions, for t/V small enough, and U + < 0, 
the model has checkerboard long range order of the electron density by following the proof 
of 1^ for the Heisenberg antiferromagnet. Here again, checkerboard long range order is 
expected to occur also in D = 2, but the proof given here does not directly apply. 

e) Extended Falicov- Kimball model. The extended Falicov-Kimball model we wish to men- 
tion has the Hamiltonian 

Hi{wx}) = E t,ycLcya + UY.hx + V (2-22) 

x,y£A,a=1,i xeA x,y£A 
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+U' (n,t + n^i - 1) {2w, - 1) (2.23) 



where {w^} is a configuration of random variables with values or 1, describing the 
position of classical particles (say nuclei or fermions with a large effective mass, we refer 
to |T^ for a discussion of the physical interpretations). The usual Falicov-Kimball model 



has U = V = and only one type of electron (say the spin up electrons). The energy of 
a nuclear configuration {wx} is Xq{H{{wx})), the smallest eigenvalue of (|2.23|) in the total 



Fock space of the electrons. A theorem of Kennedy and Lieb |T^ asserts that for the usual 
Falicov-Kimball model on a bipartite lattice A = A U B union of two sublattices A and 
B, for all U' the minimum of Xq{H{{wx})) is attained for one of the two configurations 
{wx = 0,x & A,Wx = l,x & B) OT = 0,x & B,Wx = l,x E A). This is true irrespective 
of the boundary conditions or the flux configuration (provided it exists). Many more detailed 
results are known but it is only this one that we will now generalize. 

We take a D-dimensional hypercubic lattice with periodic boundary conditions in all 
directions and set the flux configuration to be equal to tt in all square plaquettes, and also 
the canonical flux through the circuits created by the periodic boundary conditions. It is 
explained in the next section why this can be done and why it is the correct choice. We set 
\txy\ = t. By performing the sequence of transformations i)-iii) the Hamiltonian is brought 
to a refection positive form. The only term we have not discussed so far is the last one in 
( p.23| ) (the one with coupling constant U'.) After the transformations i)-iii) it becomes the 
following: 

where L and R refer to the left and right parts of the lattice with respect to some reflection 
plane P. Let r{x) denote the site obtained by reflection of x through P. It is convenient to 
use the variables = - 1 and write Eo{{sx}x^l, {sx}x€r) = >^o{H{{wx})). By applying 
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Lemma p.l| one obtains 

Eoi{s,}.eL, {s. W) > ^ (^o({s.}.gL, {s'AxeR) + Eo{{s'^},eL, {s.Wk)) (2.25) 

where for x G i? = —Sr(x) and for x G L = — By repeated reflections (across all 
reflection planes in all D directions) one concludes that the minimum energy is attained in 
the checkerboard configurations of the variables Sx (or, equivalently, of the Wx)- Note that 
this result holds for all U, U', and for all > 0. 



3 Examples and Discussions 

In this section we comment and illustrate Theorem |1.4| by various examples. First let us 
consider several planar graphs. 

1. Planar graphs. The most basic case is that of a square lattice with periodic boundary 
conditions in one direction and an even number of sites in that direction, say the horizontal 
one. Thus we have a cylinder (which can be embedded in the plane). A basic set of circuits C 
is constituted by the square plaquettes of length n = 4, and one big circle along a basis of the 
cylinder. We emphasize that if one takes only the square plaquettes then the flux through 
circuits that wind around the cylinder is not uniquely determined, and thus the set of square 
plaquettes alone is not a basic set of circuits. Obviously one can find reflection planes that 
satisfy our assumptions: these are the vertical planes that cut the cylinder in two equal 
halves. Furthermore in order to have \txy\ invariant under reflection across these planes 
we must require that {|ta;y|,< xy > horizontal} has period 2 in the horizontal direction, 
and {|txy|)< xy > vertical} is translation invariant in the horizontal direction. There is 
no constraint for \txy\ along the vertical direction. A canonical flux conflguration can be 
described by putting a flux tt through each square plaquette and 7r(A^ — 1) through the basis 
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Figure 3: A graph with non uniform optimal flux. Periodic boundary condi- 
tions are assumed in the horizontal direction. 



of length 2N of the cylinder. Theorem |LJ states that this flux configuration minimizes the 
ground state energy. 

As Lieb points out [|T], one can obtain the optimal flux on other graphs simply by "eras- 
ing", that is, letting ^ on some edges in a way that preserves the assumptions. For 
example one can get the hexagonal lattice with periodic boundary conditions where the flux 
through each hexagon is and 7r(A^ — 1) through the basis. Many planar graphs that cannot 
be obtained by this procedure satisfy our assumptions however. An example is given in 
Figure 3. 

In all these examples one can also take periodic boundary conditions in both the vertical 
and horizontal directions. This wraps the graph on a torus, thus it is not planar any more. A 
new circuit has to be added to C, namely a circuit winding around the torus in the "vertical" 
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direction. The reason for this will become clear in the next paragraph. 

2. Basic sets of circuits. Before discussing more general examples it is useful to indicate a 
way of checking that a set C is a basic set of circuits. We describe a sufficient condition. We 
start by representing any oriented circuit 7 as a sum over all edges Cj of the graph 

7 = E^i^j 
j 

The edges have a fixed reference orientation and the ej are equal to if cj does not belong 
to 7, +1 {resp. — 1) if Cj occurs in 7 with the same (respectively, opposite) orientation than 
7. We require that any circuit 7 can be decomposed as 

7 = Z!«i7i (3.1) 

i 

with integer and 7^ G C. We call such sets C generating. 

If the fiux configuration is specified for all 7i G C then using (|3.1| ) we can compute the fiux 
through 7, namely <l>^ = J2i ca^i- Since the fiux configuration through 7.j G C corresponds 
to a set of phases, different decompositions of 7 lead to the same flux. Once the flux is 
determined for all circuits, it follows from Lemma 2.1 in that C is a basic set of circuits. 

The property that C is generating can be expressed as a simple topological property of 
the surface (two-dimensional complex) consisting of the set of vertices A, the edges in F, and 
the set of triangles obtained by triangulation of all the circuits 7 G C The set C is generating 
if and only if the flrst homology group over the integers of this surface is trivial (see e.g. 
||15||). If one views the complex as a continuous two dimensional manifold this corresponds 
to the property that any closed curve can be contracted to a point. 

3. Non planar examples. For non planar graphs it is not obvious that there exists phases 
which correspond to the canonical flux. Let us consider, e.g., a single Z)-dimensional hyper- 
cube. We show that in general for a given configuration of fiuxes through the two dimensional 
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squares, one cannot find corresponding phases for tlie t^y. Tlie number of /c-dimensional sub- 
cubes is equal to 2^^^D\/{D — k)\k\. Indeed, a k dimensional subcube is determined by the 
set of points (xi, ...,xd) with < < 1, < Sj^. < 1, and xj = or 1 for j ^ ii-.-ik- 
So we have D\/{D — k)\k\ choices for ii...ik and 2^~^ choices for the Xj's. Thus the number 
of flux variables through squares is 2^^^D{D — 1), and the number of phases on the edges 
is 2^^^D. In general one will have to solve a system of equations which is overdetermined 
if 2^^^D < 2^^^D{D — 1), i.e. D > 5. However for the canonical flux configuration there 
always exist a solution of this system of equations. In fact our proof of the flux phase 
conjecture constructs such a solution, for any graph satisfying the assumptions of Theorem 

In particular the hypercubic lattice falls into our class of graphs. In order to satisfy the 
assumptions we have to take periodic boundary conditions in D — 1 or D directions. A 
generating set of circuits is constituted by all the square plaquettes and D — 1 circuits, that 
are the D — 1 coordinate axis in the periodic directions. The canonical flux configuration 
is unique and equals vr for each plaquette and 7r(A^j — 1) through the D — 1 circuits in 
the periodic directions of lengths 2Ni, i = 1,...,D — 1. One can of course imagine many 
non-planar graphs satisfying the assumptions A1-A2. 
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